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ABSTRACT: This work completes the construction of purely alge- 
braic version of the theory of non- linear quantum chemistry methods. It is 
shown that at the heart of these methods there lie certain algebras close in 
their definition to the well-known Clifford algebra but quite different in their 
properties. The most important for quantum chemistry property of these 
algebras is the following : for a fixed number of electrons the corresponding 
sector of the Fock space becomes a commutative algebra and its ideals are de- 
termined by the order of excitations from the Hartree-Fock reference state. 
Quotients of this algebra can also be endowed with commutative algebra 
structures and quotient Schrodinger equations are exactly the couple cluster 
type equations. Possible computer implementation of multiplication in the 
aforementioned algebras is described. Quality of different polynomial ap- 
proximations of configuration interaction wave functions is illustrated with 
concrete examples. Embedding of algebras of infinitely separated subsys- 
tems in algebra of the united system is discussed. 

Key ■words: Fock space; commutative algebras; non-linear wave function 
based methods 

Introduction 

At present stage of development the non-linear methods of quantum 
chemistry take on spectial significance in investigation of molecular systems 
when electronic correlation effects should be accounted on rigorous and de- 
tailed level. Among such methods non- variational version of Coupled Cluster 
(CC) method is most popular now. Its theoretic backgrounds were formu- 
lated in the seventies of the last century p] [2] and gained development in 
further numerous publications (see [3j-[10] and references therein). 

There were attempts of purely algebraic approach to non-linear theories 
of quantum chemistry (see, e.g., [H [6]) but with moderate success. And 
only after pubication of papers [IH [T2] it became clear that at the heart of 
non-linear methods of CC type there lies a very special but purely algebraic 
structure close in its definition to the well-known Clifford algebra (see, e. g., 
|151 116] ) but essentially different in its properties. And it became possible 
to give a uniform description of non-linear methods of quantum chemistry 
on modern algebraic level and compare their quality. 

In present paper we give a final reasonably rigorous report on algebraic 
backgrounds of practically all existing non-linear methods of quantum chem- 
istry and on concrete examples compare quality of different polynomial ap- 
proximations of many electron wave functions. 



Formal Fock Space 



In quantum theory notion of space of occupation numbers is widely used 
and became habitual both to physicists and to quantum chemists (see, e.g., 
|171ll8j). Recall that if the number of molecular spin orbitals (MSO) is finite 
and if they are labelled by indices from the index set A^ = {1, 2, . . . , n} then 
the corresponding space of occupation numbers is generated by bit vectors 
(rii, 712, . . .) of length n with nj = 0(1). These bit vectors constitute a basis of 
2" - dimensional vector space and their interpretation is very simple: each bit 
vector corresponds to normalized Slater determinant built on MSOs having 
non-zero occupancies in the bit vector under consideration. Support of bit 
vector is defined as a subset of the index set N, involving indices (positions) 
where occupation number is equal to 1. It is clear that each bit vector is 
determined by its support and vice versa. It is possible therefore to replace 
bit vectors by their supports and consider linear combinations of all subsets 
of the index set A^. 

Definition 1. Vector space of formal linear combinations of all subsets of 
the index set N with coefficients from, the ground field K will be called the 
formal Fock space over K associated with the index set N and will be denoted 
by the symbol 'Sn. 

Of course, in this work fields of real (K = M) and complex (K = C) 
numbers will be of actual interest. 

For each subset R C N symbol cr will stand for a basis vector of the 
formal Fock space. Of course, cr is just a notation for a subset of the index 
set A^, playing, as used here, a role of a basis vector. 

Subsets with the same number of indices generate subspace (sector) of 
the formal Fock space and 

n 

^N = ^rj, (1) 

p=0 

where 

(P) 
5^ = Ken (2) 

RCN 

With each index from the MSO index set creation and annihilation opera- 
tors can be associated and formal Hamiltonians can be defined in terms of 



these creation-annihilation operators and 1 - and 2 - electron integrals as 
parameters. 

Definition of formal Fock space immediately suggests that one can try to 
use standard set-theoretic operations (union, intersection, difference, sym- 
metric difference) to introduce in this space new laws of composition in ad- 
dition to the already existing (addition and multiplication by scalars from 
K). It should be emphasized that such an approach is widely used in mod- 
ern mathematics. Namely, first a set with some algebraic operations on it is 
taken, then vector space of formal linear combinations of elements of this set 
is constructed and initial algebraic operations are continued by linearity to 
the space of formal linear combinations. The standard example is a formal 
linear hull M.[G] of multiplicative group G. Putting 

C^O'g(^g)C^heh) = ^ aghhCgh (3) 

geG heG g,heG 

we come to the well-known and very important object called group algebra 
of (finite) group G over IK (see, e.g., [13j). 

We start with a law of composition (multiplication) on ^Af which is well- 
known for about 120 years and was first introduced by Clifford [14]. Its 
modern interpretation may be found in \15\ 116] . This multiplication is a 
combination of set-theoretic union of MSO indices and contraction of a cer- 
tain part of these indices with the aid of some fixed symmetric bilinear form 
g on '1-electron' space 5]v It is convenient to assume that basis constituted 
by the vectors e^ = era is orthogonal with respect this form. Then the 
aforementioned law of composition is defined as jl51 [16] 

Cij.es = JJ (r, s) JJ g{ei,ei)eRAs (4) 

rGR ieRns 

where RAS = RUS-RnS and 

(r..)4 \ " '•^' (5) 

1—1 if r > s 

and where it is assumed that 'empty' products are equal to 1. In terms 
of introduced in our previous publications set-theoretic operation Ak = 
AkeKi^ik] (see, e.g., [11]) the sign prefactor may be written as 

J](r,s) = (-l)l(^-^)nA(H-s)l (6) 

rGfl 

ses 



This presentation of sign prefactor is very convenient both for analytic and 
numerical purposes. For example, proof of associativity in Clifford-type 
algebras becomes practically automatic. 

Directly from definition it follows that 60 is the identity element of the 
multiplication Q . It is easy to show that this law of composition is associa- 
tive, non-commutative, and endows the formal Fock space with the Clifford 
algebra structure. Classic relation 

x.y + y.x = 2g{x,y)e0, (7) 

valid for any elements x,y G ^]^, is readily obtained from Eq.Q. 

If g is zero form then the Clifford algebra is just a Grassmann algebra 
of '1-electron' space ^\f. 

Another type of multiplications on the formal Fock space was introduced 
in our previous publications [111 [T2] . Recall briefly most important deflni- 
tions. 

For any fixed subset R C N ('reference subset') basis vectors of the 
formal Fock space may be written as ej{R) = cr-jiji where J C R, I C 
N — R. Law of composition of these basis vectors ('star product') is defined 
as (see [llj) 



e'j{R)*e';{R) 



n ir,s) 



relUJ 

sei'uj' 



e%^j, {R) if (/ U J) n (/' U J') = 0, 



(8) 







if (/ u J) n (/' u J') / 

Relations dS]) constitute the multiplication table of an algebra and it is 
easy to ascertain that this algebra is associative, skew-commutative, and 
that 60 (i?) = cr is its identity element. Skew-commutativity is readily 
follows from the equality 

e^AR) * e^XR) = (-l)l^"^l-l^'^^'leS:(i?) * e^j{R) (9) 

In the finite-dimensional case we have 2"" ways to select reference subset 
and, consequently, 2" different multiplications on the formal Fock space. Of 
course, algebras, corresponding to reference subsets with the same number 
of electrons are isomorphic. 

The simplest reference subset is i? = and the corresponding basis 
vectors are 60(0) = 6/. We have the following relations 



6^0(0)* 6^0(0) 



/- 


" 


< 


n(r,s) 

_Bei' 


h 






6^u/(0) if/n/' = 0, 
if / n /' / 



(10) 



which constitute the multiphcation table of the Grassmann algebra. In 
this particular case the star product is just the classic Grassmann exterior 
product. 

If the reference subset R belongs to the p-electron sector of the formal 
Fock space {\R\ = p > 0) and [/| = |Jl = r < min{p, n — p} then, as 
immediately follows from Eq.Q, the star product endows this sector with 
the commutative algebra structure. As will be seen later, this structure 
is particularly important for quantum chemistry. Unless otherwise stated, 
from here on only p-electron sector 5^^ of the formal Fock space will be under 
consideration. For fixed p-element subset R the commutative p-electron 
algebra with the identity eji will be denoted A^(i?). Its basis elements are 
ej{R) with \J\ = \I\ = I where I is the so-called excitation order or excitation 
level {I = 0,1, . . . ,m.m{p,n — p}). 

Algebra A^(R) admits a decomposition into a direct sum of subspaces 

A^(^) = 0v5v(i2) (11) 

where 

(0 

V5v(^)= ^e'j{R) (12) 

JCR 
ICN-R 

and where Yj^{R) = {0} for I < and / > min{p, n — p}. For fixed 1 < / < 
min{p, n — p} this subspace can be called 'subspace of all excitation of order 
/ from the reference subset R\ From the definition of the star product it 
readily follows that 

V^(i?)*v5v(i?)cV^+'(i?) (13) 

which means that A^{R) is a "L-graded algebra. 
Subspace 

rmn{p,n—p\ 

iUR)= v^(i?) (14) 

l=r+l 

is an ideal of A^{R) for any r > 0. This means that any quotient space 
A^{R)/Tj^{R) is also a commutative algebra with the law of composition 

[^1 * ^2]r- = [*l]r. * [*2]r (15) 

where [^]r = ^ + I^(i?). By an abuse of language star product in the 
quotient algebra A^(i?)/I^(i?) will be called 'star product modulo r'. 



It is clear that I^(-R) is a maximal nilpotent ideal of A^{R). Maximality 
means that it is not contained in any ideal of the algebra under consideration 
and nilpotency means that [l%{R)Y = 1%{R) * . . . * 1%{R) = {0} for r > 
min{p, n — p}. 

A^(R) is an abstract model of a typical full Configuration Interaction 
(CI) space (endowed with the additional multiplication) devoid of any orbital 
specificity. Vectors from the vector space 

r 

W5v(ii) = 0V^(i?) (16) 

1=1 

will be called amplitudes. Note that in contrast to the standard CC type 
methods our amplitudes are just wave functions from the same p - electron 
sector of the Fock space. 

If J-^ is a concrete 1-electron space and ip = {ipi,ip2, . . . ,ipn) is some its 
orthonormal n-frame of MSOs, then the substitution mapping 

s^-.CR^ Vii A ^i2 A ... A Vip, (17) 

where R = {21,^2 .. . ,«p}, ii < 12 < ■ ■ ■ < ip establishes (basis-dependent) 
isomorphism between the formal Fock space '^m and concrete Fock space 

n p 
p=0 

Substitution mapping can also be used to transfer the algebra structures 
defined by Eq.([8|) from the formal Fock space to any its concrete realization. 
The inverse mapping si (also basis dependent) strips vectors of concrete 
Fock space of their orbital specificity. 

It pertinent to emphasize that for any concrete realization of the Fock 
space star product of arbitrary p - electron wave functions is (antisymmetric) 
p - electron wave function. 

The formal Fock space can be endowed with Euclidean (K = M) or 
Hermitian (IK = C) structure in a simple and natural way: 

, , ^ (1 [fR = R' 

If a concrete 1-electron Fock space is a vector space with the inner product 
and f/' is some its orthonormal (with respect to this inner product) n-frame 
then the substitution mapping is obviously isometric. 



The group of orthogonal (or unitary) transformations of 1-electron Fock 
space T^ wiU be denoted Q(T^). If ro/j isp-plane (p-dimensional subspace) of 
1-electron Fock space spanned by MSOs {"01 }> i ^ R and wj^ is its orthogonal 
complement then for any block-diagonal transformation U G Q{wji)xQ{wj^) 

s^V'^iR)) = s^f/(V^(i?)) (20) 

where I = 0, 1, ... , min{p, n — p}. In the case / = this equality is equiva- 
lent to the well-known property of the determinant states: each determinant 
state is defined up to arbitrary orthogonal (unitary) transformation of occu- 
pied (virtual) MSOs. Therefore, Eq. (|20p may be considered as generalization 
of this simple property of determinants. It is to be noted that, in general, 
s^(r) 7^ S0'(t) for concrete non-zero vector r G Vjy(-R)- 



Polynomial Functions With Algebra A^(i?) As Their 
Range 



Many methods of quantum chemistry may be uniformly described in 
terms of polynomial functions Pa : W^(i?) — > A^(i?) of the form 

min{p,n—p} 

Pa{r)= Y. ^^^' (21) 

/i=0 

where Cj € IK and t^ = t * ■ ■ ■ *t . Since W^(i?) is a subspace of the 

k times 

nilpotent ideal 1^{R), all polynomial functions of the type of Eq. (j21l) have 
degrees not greater than min{p, n — p}. 

Using polynomial functions it is easy to characterize invertible elements 
of algebra A^{R). 

Propsition. Wave function ^ is invertible element of algebra A^{R) if 
and only if ao = (^|e|(i?)) / 0. 

Proof. Let us put r = 'I' — aQeQ{R). By direct calculation it is easy to 
ascertain that 

mm{p,n-p} ^ 

^-'= E (-1)"4+T (22) 

At=0 "O 

Corollary. Amplitude vectors are not invertible. 



Polynomial functions PaiT) such that (1) ao 7^ 0, and (2) PaiT) is injec- 
tive, are most important for applications in quantum chemistry. The first 
condition means that for any amplitude vector r the state vector Pa{T) has 
non-zero projection on the reference basis vector e0{R), and, consequently, 
is invertible. Without loss of generality we can put oq = 1- The second 
condition guarantees that the mapping r -^ PaiT) is a parametrization of 
the surface S'^{R) = Pa (W^(i?)) in the algebra A^(i?) of p-electron states. 
Dimension of this surface is equal to the dimension of the amplitude space 

w^R)- 

Propsition. Inverse of polynomial function PaiT) ■ W^(i?) —?■ S'^iR) 
is also a polynomial function and it exists if and only if ai 7^ 0. 

Proof. Follows directly from the recurrence formulas for the coefficients 
of the inverse polynomial P~^ iaQCQiR) + t) = ^ (— 1)^^^6^t^: 

&M = ^^E(-1)'"'^- E «n---«^., (23a) 

hi = —. (23b) 

Note that invertibility of PaiT) is equivalent to the condition that the 
polynomial derivative 

(min{p,?i— p} \ min{p,n— p} 

E «M^M = E ^'^^^'~' (24) 

is non-zero for any amplitude vector r as could be expected. 

For any reference subset R and any polynomial function PaiT) with 
non-zero oq and oi a triple cr = (5]y(i?), [Pa]""*^, W^(i?)) is a chart of a 
geometric object which can be called polynomial manifold. Note that any 
point of the surface S^iR) is invertible in algebra A^(i?). If wave function 
^ belongs to the surface 5]y(i?) then ^~^ belongs to the surface which can 
be called 'surface of inverse elements of 5]y(i?)' and, by a certain abuse of 
language, will be called surface inverse to SJ^iR). Thus, each polynomial 
surface (manifold) arises in pair with its inverse. 

It is pertinent to note that for a fixed polynomial parametrization all 
arising surfaces form a chain 

ShiR) c sUR) c. . . c 5^^"^^'"-^^(i?) (25) 



and, consequently, any surface S'^{R) can be treated as a subsurface of 
S^{R) for any r' > r. As a result, it is possible to take some reference tq 
corresponding to a fixed linear combination of determinants and involving 
certain excitations of arbitrary order from the reference subset R. Consid- 
ering PaiTo) as the origin, it is possible to construct parametrizations of the 
type r i-)- Pa{To)*Pa{T) where r include all excitations up to preferable maxi- 
mal order r. If tq = 0, we come to the standard polynomial parametrization. 
Such an approach presupposes that tq is sparse to be kept in fast memory. 
Among infinitely many polynomial functions at present only three types 
are of actual use in quantum chemistry. The first type is certainly the affine 
function 

Pci:t^ e^{R) + r, P^} : e%{R) + t ^ t (26) 

which corresponds to the CI approach accounting all excitations up to order 
r from the reference index set R. The corresponding surface SC^{R) is just 



the affine plane eS (R) + W^ (R) . The inverse surface is 



min{p,n— p} 

r^ E (-l)V' (27) 

fe=o 

The second type is the exponential function 

min{p,ra— p} min{p,n— p} 

exp:r^ ^ ^, exp'^ : e0(i?) + x ^ ^ (_i)M-i^(28) 

which appears in coupled cluster methods. Exponential surface SE^j^{R) 
coincides with its inverse and is, in fact, a multiplicative Abelian group. 
And the last, third type is 



Qa ■ T 1-^ eQ{R) + T + aT'^, q^^ : eQ{R) +x ^ x - ax^ -\ x^ H 

(29) 
corresponding to the quadratic CI (QCI) method p^. With such a parametriza- 
tion calculation of polynomial value can be done very fast but to reach rea- 
sonable precision in approximation of wave function excitation level should 
be sufficiently high. 

Parametrization closely related to the resolvent mapping in algebra A^{R) 
can also be used 

min{p,n— p} min{p,n— p} 

Q:t^ Y, r\ Q-^:e%{R) + x^ ^ (-l)'^-V' (30) 

i=o ^l=l 

10 



It is easy to see that surfaces SQ^j^{R) and SC'^{R) are mutually inverse. 

If r = mm{p, n — p} then any injective polynomial function may be used 
to parametrize all p-electron states that have non-zero projection on the 
reference basis vector e0{R). This means that from the variational view- 
point all the aforementioned polynomial parametrizations are equivalent. If 
r < ram{p, n — p} then different polynomial functions parametrize different 
low dimensional surfaces in the space of all p-electron states and correspond, 
in general, to different variational methods. In this case the choice of ap- 
propriate polynomial parametrization is of primary importance. 



Variational Methods on Polynomial Surfaces 



In this section, unless otherwise indicated, it will be supposed that the 
reference index set R is fixed, and symbol R in all expressions will be su- 
pressed. 

For arbitrary polynomial function Pa{T) the electronic energy expression 
may be written in two equivalent forms: 



{Pa{s4T))\H\Pa{s4T))) 

\Pa{r) 



■^'^('^) = IIP ,,M|2 (31a) 



{Pa{T)\H^\Pa{T)) 

where ip is arbitrary n-frame of J-^, H is the standard p-electron Hamil- 
tonian, and H^ is its formal parametric analogue expressed via formal 
creation-annihilation operators 

n ^ n 

i,j=l i,j,k,l=l 

Remind that s^ is isometry and, consequently, |li^a('5v (''")) IP = ll-fa(''")lP- 
Unless otherwise indicated, the energy expression (I31b[) will be used. 
To get energy stationary conditions, it is necessary to find out interre- 
lation between star product and derivation. For K = M electronic energy 
is analytic function of amplitudes whereas for K = C it is not. However, 
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realification of the complex amplitude space improves the situation. Let us 
put 

r (0 
'^ = Y. Y. {x'j + w'j)e'j = x + \y (33) 



1=1 JCR 
ICN-R 



Then 



and 



^Pa{r) = P'At) * eS, -^Pa{r) = iP'air) * e'j (34) 



^(^) = ||p^^^^p Re(P^(r) * e5|F^ - Ea{T)I\Pa{T)) (35a) 

9y5 lKa('r)||^ 

where the expression for polynomial derivative P'a{T) is given by Eq. (l24p . 

Now we can write down stationary conditions for electronic energy as a 
function on the surface 5]^: 

(P^(r) * e'j\H.^ - Ea{T)i\Pa{T)) = (36) 

where J C R, I C N - R, \J\ = \I\ = 1,2, . . . ,r. 

Wave functions Pa{T)*ej and iP^(r)*ej are linearly independent forming 
a basis of the realified tangent space Tp mS^ to the surface S^ an the point 
Pair). 

In the particular case Pa{T) = e^ + t system (f36l) is equivalent to the 
eigenvalue problem for the projection of the electronic Hamiltonian on the 
subspace spanned by the basis vectors Cj with < |/| = | J| < r (CI method 
accounting all excitations up to order r from the reference determinant). 
In general case the system (136p can be reduced to eigenvalue problems for 
electronic Hamiltonian projection on r - dependent subspace spanned by 
vector Pai^) and tangent vectors Pai^) * ej. Since these vectors are not 
orthogonal, non-trivial Gram (overlap) matrix arises and the aforementioned 
eigenvalue problem is actually a generalized one. 

At this stage we can see advantages of the exponential parametriza- 
tion. First, calculation of the polynomial derivative is not required because 
exp'(r) * ej = exp(T) * Cj. Second, algorithm for r update is very simple 
and efficient and can be described as follows. 

(i) For a fixed amplitude vector r generalized eigenvalue problem for 
Hamiltonian projection on subspace spanned by vectors expr * Cj with < 
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|</| = l-^l < ^ is solved to give the ground state vector 

r (0 

^(r) = coexp(r) + J]] ^ c5exp(T)*e5 (37) 

1=1 JCR 

ICN-R 

(ii) After dividing by cq wave function is rewritten in the form 

^(r) =exp(T + Ar) (38) 

where 



(0 



\ 



Ar = lo9\e^ + J2 E ^Wj 



= 1 JCR 

ICN-R 



(39) 



/ 



(iii) If 1 1 At 1 1 is still greater than some threshold value, we put t -i^ t + At 
and return to step (i). 

At each step of this algorithm energy value never rises. Note that at 
step (ii) we used the fact that exponent is an isomorphism of additive group 
of amplitudes onto multiplicative group SE^. For arbitrary polynomial 
parametrization the analogous algorithm is more complicated. 

In the same simple manner it is possible to derive linearized equations for 
system evolution on arbitrary polynomial surface in the vicinity of electronic 
energy stationary point. To this end we should first get formula for electronic 
energy second derivatives assuming that the ground field is C Realification 
of the amplitude space and calculation of the second energy derivatives with 
respect to real variables x, y give 



D^Eaix,y) 

W-faV )\\- 

_Tm r/'TT J- A ^ _ p rn J- t;M \ 

(40) 

^ ±111 ^y±± - ^j - ±:j yKj - u;j iLC iyi^± - l^ j - jj. yKj - u j i j 

where 



Re [(H + A) - S (G + S)] -Im [(H + A) - E (G + S)] 
Im [(H - A) - ^ (G - S)] Re [(H - A) - E (G - S)] 



Hj/,j'r = (P^(t) * e'j\H\P',{T) * eS',) (41a) 

Gj7,j'/' = kKij) * e'j\K{r) * eS',) (41b) 

Aji,.rv = {Pa{T)\H\P':{T) * e'j * eS',) (41c) 

Sj/,j'/' = {Pa{T)\P':{T) * e'j * eS',) (41d) 

Here H is a Hermitean matrix representing one of two identical blocks of 
projection of the electronic Hamiltonian on the realified tangent space to 

13 



the surface S^ at a point Pa{T), G is a Hermitean Gram (overlap) matrix 
of, in general non-orthogonal, basis {P'(t) * ej}, A and S are symmetric 
matrices arising due to non-zero curvature of the aforementioned surface. 

Using standard technique based on the theory of Hamiltonian equations 
on symplectic manifolds (see, e.g. [22]- [25])' i* is easy to derive the following 
linearized evolution equation in variables r, f (Schrodinger form) 

H-K(ro)G A-S J , , ( ^^2) 



" \\Pa{roW V -(^ + S) -(H - K(ro)G) 

where tq is electronic energy stationary point, matrices H, G, A and S 
are calculated at this point, and the amplitude vector r in this equation 
is just the replacement vector from tq. In quantum chemistry such an ap- 
proach to calculation of transition energies is called time-dependent (TD). 
Thus, system (j42p embraces all versions of TD methods based on polynomial 
parametrizations, including all versions of TDCC and TDQCI methods. 

In simple case of parametrization Pa (r) = eS + r Gram matrix reduces 
to the identity matrix, matrices A and S just vanish. And as could be 
expected, TDCI method, due to zero curvature of the surface 5C]y, gives 
nothing new in compare with the linear CI approach. 

In concluding this section it is pertinent to give a very brief outline of 
possible generalization of variational problem under discussion. If MOs are 
to be varied, then the electronic energy domain may be elegantly described 
as a (locally trivial) vector bundle with Grassmann manifold over 1-electron 
sector of Fock space as its base and formal amplitude vector space as a typical 
fibre. Such an approach embraces Hartree-Fock (HF) case (zero typical 
fiber) and Multi-Configurational Self-Consistent Field (MCSCF) methods 
along with already discussed ones, and all relevant TD methods. 

At present stage of development fully variational non-linear methods are 
considered as unfeasible. The reason is in problems with Hamiltonian matrix 
element evaluation. Even if the current amplitude vector r is of reasonably 
small length to be kept in fast memory, full CI expansion of Pa{T) required 
for matrix element evaluation may be of huge size and level of complexity 
of matrix element evaluation is close to that in the FCI method. One of the 
approaches which may help in avoiding this bottleneck is the following. It is 
possible to get complicated but closed analytic expression for vector HPa^r) 
without preliminary FCI expansion of Pa{T). Arising on this way overlap- 
type matrix elements can be handled recurrently also without explicit CI 
expansions. However, this is only a hint for possible future developments. 
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Polynomial Functions With Quotient Algebras As 
Their Range 

For each fixed excitation level r and each polynomial Pa '■ WJ^{R) — )• 
A%{R) we can consider the function P^ : W^{R) -^ A^(i?)/I^(i?) defined 
as 

P^-.T^ [Pa{T)]r (43) 

Using star multiplication in quotient algebra, we can write 

r r 

^:w = e«a«[^^]'- = E«mH^ (44) 

Note that for calculation of polynomial value in quotient algebra only the 
first r + 1 coefficients a^ are required. If ai ^ then the inverse polynomial 
[P^']~^ exists and its coefficients bi,b2, ■ ■ ■ ,br can be calculated using general 
Eq.(l23]) (60 = by definition). 

Since Fock space is equipped with the inner product (see Eq. (fT9|) ). we 
can identify quotient space A^(i?)/I^(i?) with the orthogonal complement 
to I^(-R) which is just the subspace Ke^{R) WJ^{R). 

Strictly speaking, for a fixed polynomial parametrization PaiT) the Schrodinger 
equation 

HPair) = EPair) (45) 

has a solution if and only if 

(i) the exact eigenfunction of the electronic Hamiltonian H has non-zero 
projection on the chosen reference subset R; 

(ii) the excitation level is sufficiently high (ideally r = min{p, n — p})- 

In other words, the eigenvector should belong to the polynomial surface 

5]y. In practice, however, relatively low excitations levels are used. This 

means that we can hardly hope to find amplitude vector, satisfying equation 

(jl5]l . But it is possible to consider quotient Schrodinger equation in quotient 

[HPa{r)]r=E:{r)[PaiT)]r (46) 

which can be recast in the following more habitual form 

{e'j\HPair))=E:{r){e'j\Pa{T)) (47) 

{J CR,I CN -R,0<\J\ = \I\<r) 

The number of scalar equation in this system is dim[PV^(i?)] + 1, that is 
exactly the number of unknowns (amplitude vector coordinates) plus 1 . One 

15 



additional equation corresponds to J = / = and is used to evaluate 
'energy' value for a fixed r: 

Kir) = ^ 0' "^ '^ (48) 

With such an approach 'energy' E"^ depends both on the paranietrization 
(polynomial coefficients a) and excitation level (r), not being variational. It 
tends to the exact value when r -^ min{p, n — p} independently of a. 

Iterative solution of linearized system (I47p is normally used in CC type 
calculations (see, e.g., [26]). One can try as well to solve iteratively the 
initial vector equation (I46p rewriting it in the form 



l/ [ffP,(T)]A , [HPa{T)]r 



[HPa{r)]r ' 
Kir) 



+ ■■■ (49) 



In concluding this section it is pertinent to note that, when solving the 
aforementioned equations, the main problem arising is the calculation of 
vector [HPaiT)]r required at each iteration. To perform such a calculation 
CI expansion of the polynomial function PaiT) up to excitation level r + 2 
is required. Using star product, we can instead derive formulas for direct 
action of creation-annihilation operator products on polynomial functions. 



Computer Implementation Of Star Product In Quo- 
tient Algebras 



When using polynomial parametrization, one of the crucial points is the 
calculation of star product of two wave functions modulo r. Since almost 
all calculations in quantum chemistry are performed in MO representation, 
from the very beginning we separate a - and /3 - indices and replace MSO 
index sets by pairs of MO index sets. In particular, MSO index set A^ 
will be replaced by a pair {M,M) of a - and /3 - index sets, |M[ = m, 
reference subset R will become a pair (i?a, Rp), and basis elements of algebra 

A,j^'j^^{Ra,R(s) will be efjj {RajRp). In what follows, unless otherwise 
indicated, it will be supposed that the reference pair of index sets {Ra,R^) 
is fixed, and in all expressions dependence on {Ra.,Rp) will be suppressed. 
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Let us present classes of wave functions modulo r to be * - multiplied in 
the standard form 

fc=0 ka JaCRa Jj^CRp 

laCM-Ra /^CM-flo 

r (ka) (kis) 

['^]^ = i^ + EE E E 4:4yti: (50b) 

A; = ka JaCRa Jfi^R^ 

laCM-Rc, /^cM-flo 

where ka + kp = k and max{0, k—pp, k — m+pp} < ka < minjpQ, m—pa, k}. 
One of possible expressions for star product of two wave functions is 

r (fca) (fc/y) 

E 

(51) 



fc=0 



ka JaCRa Jp^Rjl fc'=0 fc^ J'^CJa J'^-ZJ^ 



where k'a + k'a = k' , max{0, k' — kjs} < k'^ < min{/c', ka}, and 

e= Y, [|(j;u/;)nA(,^u/.)l + ^;] (52) 

cr=a,(5 

Summation limits over ka and k'a require probably some comments. We 
have inequalities 

< /cq, < Pq 

< ka < m — Pa 

< k - ka < Pi3 

0<k — ka<m— pi^ 

and, consequently, ka > max{0, /c— p^, fe — m+p^}, ka < Tain{k,Pa,m—Pa}- 
As for k', we have 



"q; 



< fe^, < fe« 
< k' - kL < ka 
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which means that max{0, k' — k/^} < k'^ < min{fc', ka}- Simple part of star 
product in Eg . ([51]) corresponds to k' — xj, rv auLi id -^0 0-- 
The second version of formula for star product is 



0, k and is x^'®$ + *2/g'g. 



r (ka) 



EE 



{K) 

Jq, C-Rq: — J a 
ILCM-Ra-Ia 



fc=0 fcc 



ik'e) 






i.kp) 

E 



J en J 13 ' 



laCM-Ra Ia(ZM-Rp 



J'„CR 



JcUl'aJfiUr^ 



r' T' 



■13 -■'13 



e= Y. |(j;u/;)nA(j^u/. 



(53) 



(54) 



a=a,li 



In the above equations r - maximal excitation level used. For full (not 
quotient) product r = min{p, 2m — p}- Note that simple contributions in 
the last expression correspond to A; = and k' = 0. 

At present we are using the second expression for the star product mod- 
ulo r because in our realization the corresponding algorithm works faster 
than that based on the first expression. 

Any computer implementation of start product requires algorithm for 
generation of subsets with fixed number of elements contained in a given 
index set. For this there exist variety of methods which differ in subset or- 
derings: lexical, Gray, revolving door, etc. Apart from generation itself, very 
important also is availability of the corresponding ranking and unranking 
functions. 

After series of experiments we came to the conclusion that the revolving 
door algorithm is the best for our purpose. This algorithm and accompa- 
nying ranking and unranking functions are described in [27J. The corre- 
sponding codes are free and can be found at The Stony Brook Algorithm 
Repository [28j. 

For our current purpose unranking function is of no use, but ranking 
one is compulsory. Let us denote by rank^j J rank of J as a subset of R in 
revolving door ordering, and by rank(^^ ^^ )(Jq,, J/3) rank of a pair of subsets 
calculated as 



rank(K^_^^)(Ja, J^) = (rank/j^ J„ - 1) I , I + rank/j^ J/3 



Pa 



(55) 
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Then for given k, ka, Ja, Jp-, la-, and Ip the position of the corresponding 
element in a row vector representing wave function (as well as its amplitude) 
is determined with the aid of the following mapping: 

{k,ka,Ja,J(3,Ia,I(3) ^ oEset{k, ka) + 

(56) 

and where the array of offset indices (filled only once at the beginning of 
calculation) is 

fc =1 fe^=max|0,K — p;3,fe — ?Ti+p^| ^ '"^ 

Efpa\ fm - Pa\ fpjA fm - Pj3 
\k'J\ k' J\k')\ k' 

(57) 

Computer code was written in FORTRAN 95 and was heavily tested. Ex- 
ample of time characteristics of current computer implementation of star 
product modulo r are given in Table [H 
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Excitation 


Dimension of 


Time 


level 


the amplitude 


(sec) 


r 


space 




3 


18818 


0.02 


4 


98693 


0.35 


5 


294965 


2.75 


6 


558809 


13.04 


7 


755081 


38.49 


8 


834956 


74.62 


9 


851956 


103.10 


10 


853702 


115.27 


11 


853774 


117.95 


12 


853775 


118.23 



Table 1: Timing of star product modulo r calculations (Intel Fortran, CPU 
Q9550@2.83GHz, sequential mode). 
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Comparison of Different Polynomial 
Parametrizations On Concrete Examples 



Let us suppose that for a molecule under consideration orthonormal 
MO set is generated at the HF level of theory, active space is chosen, and 
FCI calculation of few low-lying electronic states is perforined in this active 
space. Let ^fci be one of the calculated wave functions. Our nearest aim 
is to compare this exact wave function with its approximation by different 
polynomial functions for different excitation levels. To this end we first 
use GAMESS US program p9] (ALDET route) to generate coefficients of 
FCI expansion over determinant basis set. Then maximal by absolute value 
FCI coefficient is found and the corresponding determinant is taken as the 
reference one. Next step is to reorder FCI coefficients in the ordering defined 
by the mapping (j56|) (in GAMESS lexical ordering is accepted). Then FCI 
wave function is divided by the coefficient cq = Cqq and for a selected 
polynomial parametrization the exact amplitude vector is determined: 

Tr'' = p-He^% + x) (58) 

where x = ^*fc/ - e^'^. 'Exact' means that PaiT^""""^) = ^^FCi- This 
amplitude vector depends on the parametrization used, carries therefore sub- 
script a, and involves, in general, all possible excitations from the reference 
state. 

There are two questions arising. 

(i) Suppose that the exact amplitude vector is truncated to the level r, 
and full expansion of the wave function PaiT) is calculated. The question 
is: At what excitation level polynomial of exact but truncated amplitude 
becomes reasonably close to the FCI function? 

(ii) Suppose that truncated amplitude vector is optimized to make ||Pa(r)- 
^FC/|| as small as possible. The question is: In what cases optimization may 
essentially improve polynomial approximation of ^ pci'^ 

Of quite a large number of performed atomic and molecular calculations 
in different bases and with different active spaces we give here the results 
for only three molecular systems: (1) calculation of hypothetic molecule 
H2C — NH — CHF ground state which is strongly correlated; (2) calculation 
of N2O molecule ground state; (3) calculation of the first excited (triplet) 
state of boric acid molecule (its ground state is of no interest for us since it 
is practically non-correlated). 
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Polynomial Excitation Dimension of the 
Pair) level amplitude space 



\Pa{T)-^FCl\\ \\Pa{T*)-^FCl\ 



exp(r) 



1 
2 
3 
4 

5 

6 



64 

1424 

12624 

55324 

135068 

208764 



0.420632E+00 

0.153643E-01 

0.207110E-02 

0.123119E-03 

0.114390E-04 

0.375204E-06 



0.417226E+00 

0.144873E-01 

0.191069E-02 

0.113253E-03 

0.114390E-04 

0.375204E-06 



Q{r) 



1 
2 
3 
4 
5 
6 



64 

1424 

12624 

55324 

135068 

208764 



0.438430E+00 

0.550005E-01 

0.270095E-01 

0.183023E-02 

0.746502E-03 

0.118016E-04 



0.425482E+00 

0.405679E-01 

0.194293E-01 

0.135976E-02 

0.532138E-03 

0.888355E-05 



Pair) 



1 
2 
3 
4 

5 

6 



64 

1424 

12624 

55324 

135068 

208764 



0.406908E+00 

0.396744E-01 

0.249536E-01 

0.112012E-02 

0.709024E-03 

0.872110E-05 



9i(r) 



1 

2 
3 
4 
5 
6 



64 

1424 

12624 

55324 

135068 

208764 



0.420542E+00 

0.163069E-01 
0.390887E-02 
0.157044E-02 
0.917063E-03 
0.438127E-04 



0.417171E+00 

0.151996E-01 
0.359143E-02 
0.139820E-02 
0.815560E-03 
0.389069E-04 



Table 2: Comparison of different parametrizations: Ground state of molecule 
H2C - NH - CHF. 
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Polynomial Excitation Dimension of the 
Pair) level amplitude space 



\Pa{T)-^FCl\\ \\Pa{T*)-^FCl\ 



exp(r) 



1 
2 
3 
4 

5 

6 



72 

1818 

18818 

98693 

294965 

558809 



0.415883D+00 

0.110642D+00 

0.415729D-01 

0.938115D-02 

0.803279D-03 

0.450764D-04 



0.400344D+00 

0.106583D+00 

0.392551D-01 

0.887311D-02 

0.744422D-03 

0.412938D-04 



Q{r) 



1 
2 
3 
4 
5 
6 



72 

1818 

18818 

98693 

294965 

558809 



0.628926D+00 
0.347899D+00 
0.164148D+00 
0.599690D-01 
0.168746D-01 
0.191710D-02 



0.4495243D+00 

0.1994256D+00 

0.9625347D-01 

0.3734982D-01 

0.1129635D-01 

0.1073262D-02 



Pair) 



1 
2 
3 
4 

5 

6 



72 

1818 

18818 

98693 

294965 

558809 



0.369100D+00 

0.171212D+00 

0.727427D-01 

0.211776D-01 

0.435432D-02 

0.556137D-03 



9i(r) 



1 

2 
3 
4 
5 
6 



72 

1818 

18818 

98693 

294965 

558809 



0.405992E+00 

0.119679E+00 

0.822454E-01 

0.685212E-01 

0.446974E-01 

0.135399E-01 



0.394752E+00 

0.110201E+00 

0.579319E-01 

0.418136E-01 

0.281541E-01 

0.771915E-02 



Table 3: Comparison of different parametrizations: Ground state of molecule 
N2O. 
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0,030 
0,025 
0,020 
0,015 
0,010 
0,005 
0,000 




Molecule H^C-NH-CHF: Ground state 
Basis 6-31 G* 

Number of inactive IVIOs - 1 2 
Number of active IVIOs -12 
Number of correlated electrons - 8 
Dimension of FCI space - 245025 



FCI limit 



Excitation level 



Figure 1: Comparision of different polynomial approximations on example 
of the ground state of molecule H2C — NH — CHF. Curves corresponding to 
optimized amplitudes carry additional label *. 

Comparision of quality for different polynomial appriximations of FCI 
functions for the aforementioned molecules are depicted in Figs.([l|)-(l3]) and 
presented in Tables ([I])-®. From Figures and Tables it is seen that all 
polynomial parametrizations give reasonably good approximations of FCI 
functions starting from the excitation level 6. For lower excitation levels 
exponential parametrization is beyond comparision giving the best results. 
Behaviour of polynomial approximations in the course of amplitude opti- 
mization is of special interest. Such optimization may essentially improve 
wave function approximation for, say, resolvent {Q polynomial) or quadratic 

(gi polynomial) parametrizations. In the CI case optimal truncated am- 

2 
plitude vector just coincides with the truncated part of the CI expansion 

divided by cq. And it seems that the exponential parametrization inherits 
this property of truncated CI expansions. Namely, it possesses the following 
for the time being empirical but important property: 

For exponential parametrization optimal amplitude vector involving ex- 
citations up to level r is a very good approximation to the exact amplitude 
vector truncated to the same excitation level. 
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IIPaM 

0,18. 
0,17 
0,16 
0,15. 
0,14 
0,13 
0,12. 
0,11 
0,10 
0,09 ■ 
0,08 
0,07 
0,06 ■ 
0,05 
0,04 
0,03 ■ 
0,02 
0,01 
0,00. 
-0,01 




Molecule N^O: Ground state 
Basis aug-cc-pvDZ 

Number of inactive MOs - 5 
Number of active MOs - 12 
Number of correlated electrons - 12 
Dimension of FCI space - 853776 



FCI limit 



10 12 

Excitation levels 



Figure 2: Comparision of different polynomial approximations on example 
of the ground state of molecule N2O. Curves corresponding to optimized 
amplitudes carry additional label *. 

It is to be noted that quadratic parametrization gives, as a rule, ap- 
proximations of wave functions close in their quality to the exponential 
parametrization. However, in a number of cases behaviour of such sim- 
ple parametrizations may be strange and unpredictable (see Fig. ([2])). It is 
therefore necessary to keep in mind that quality of simple parametrizations 
(quadratic, qubic, etc.) may be system dependent. 

And in concluding this section it is pertinent to mention that all parametriza- 
tions with excitation levels 1 and 2 in the amplitude vector give poor ap- 
proximations of wave functions as seen from Tables ([Z])-®. 

Wave Function of Non-Interacting Subsystems 



We confine ourselves to the case of two non-interacting subsystems. 
Let us assume that the first ni vectors e.j span MSO space of subsystem 
I, and the remainder vectors Cj span MSO space of subsystem II. The 'one- 
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0,005 



0,004 



0,003 



0,002 



0,001 



0,000 



CI 




Molecule B(0H)3:1st triplet state 
Basis aug-cc-pvDZ 

Number of inactive IVIOs - 10 
Number of active IVIOs - 13 
Number of correlated electrons - 12 
Dimension of CI space - 2944656 

FCI limit 



Excitation level 



Figure 3: Coniparision of different polynomial approximations on example 
of the 1st triplet state of boric acid molecule. 



electron' sector of formal Fock space is a direct sum 



^N — 5^Xi ® S^ 



>X2 



(59) 



where Xi = {1, 2, . . . , ni} and X2 = {rii + 1, ni + 2, . . . , n}. 

Let ^/ G 5"^ and "^u G ^^ . Wave function describing system of 
two non-interacting subsystems is just the Grassmann product */+// = 
^/ A ^// G 5^^ ^^. Is it possible to present such a wave function as a 
star product of ^/ and ^//? The answer is positive. Indeed, let us select 
reference subsets Ri C Xi and R2 C X2 for states ^7 and ^//, respectively, 
and expand these states in a standard manner 



rnin{pi,ni— pi} 



/l=0 -hCRi 

IlCXi-R-i 



mm{p2,n2—p2} 

^niR2)= E E 44(R^ 



(60a) 



(60b) 



l2=0 J2CR2 

l2CX2~R2 
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Since, by definition, ej_{Ri) = e^ ' '(0)) we can write 

ggi-JiU/i(0) A eg2-^2^^2(0) = eg^^u^(0) (61) 

where R = R1UR2, J = Ji U J2, and / = Ji U/2. The sign on the right-hand 
side of the last equality is always plus due to the chosen ordering of MSO 
indices. But due to the same ordering we also have 



e 



R-JUlfrA\ — Jl ( T)\ ^ J2 





i0) = e'j\iR)*e%iR) (62) 



Immediate consequence of this relation and multiplication law ([5]) is the 
following equality 

^i{Ri) A ^niR2) = ^liR) * ^iiiR) (63) 

In fact, in the last equality we used implicitely the embedding of g - electron 
algebra A^(S') into p - electron algebra A^(i?) {q <p). 

With words the equality ()63p can be explained in the following way. It 
is necessary to take wave functions describing non-interacting systems and 
expand them using their reference subsets and their MSO index sets. Then 
each of these functions should be embedded into p - electron sector of the 
Fock space where united reference set Ri U R2 is chosen, and then star 
product should be calculated. 

Conclusion 



This work is completing the construction of general algebraic theory of 
non-linear methods of quantum chemistry which was started in our pre- 
vious publications. Efficient computer implementation of star product for 
the most general case made it possible to compare different polynomial ap- 
proximations of many-electron wave functions for different excitations levels. 
Further possible generalizations embracing theories with MO optimization 
require essentially more complicated mathematical tools of modern differ- 
ential geometry, in particular, theory of vector bundles. On this route it is 
possible to construct theory including practically all existing wave function 
based variational methods of quantum chemistry. 

But in the frameworks of described in this work purely algebraic ap- 
proach there still exits one unsolved fundamental problem. This problems 
concerns the question of consistency between star product and the standard 
Hermitean norm on the Fock space. In more rigourous terms, it is desirable 
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to ascertain that the inequality ||^*<5|| < ||^||||$|| holds true for arbi- 
trary p - electron wave functions when the number of MSOs is sufficiently 
large. If this statement is correct, then practically all the results described 
in this work can be easily reformulated for infinite-dimensional Fock spaces 
and then the star product will become the property of quantum mechanics 
but not only of its finite dimensional model which is at present quantum 
chemistry. 
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